A mechanism for the non-Fermi-liquid behavior in CeCug a Au x 
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We propose an explanation for the recently observed non-Fermi-liquid behavior of metallic alloys 
CeCu6-a; Au^: near x = 0.1, the specific heat c is proportional to Tln(Tb/T) and the resistivity 
increases linearly with temperature T over a wide range of T. These features follow from a model in 
which three-dimensional conduction electrons are coupled to too-dimensional critical ferromagnetic 
fluctuations near the quantum critical point, x c — 0.1. This picture is motivated by the neutron 
scattering data in the ordered phase (x — 0.2) and is consistent with the observed phase diagram. 
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Fermi-liquid (FL) theory has traditionally led to an 
accurate description of the low temperature properties 
of metals. Even in the heavy-fermion compounds, where 
the bare electron parameters are renormalized by up to 
three orders of magnitude by the interaction, FL behavior 
is observed at low temperatures T with a specific heat c oc 
T, a magnetic susceptibility x ~ const, and a resistivity 
p ~ Po + AT 2 . However, in several heavy-fermion systems 
IjJ-Q pronounced deviations from FL behavior have been 
found in a number of physical properties. 

Three main theoretical scenarios have been proposed 
to explain the occurrence of the non-Fermi-liquid (NFL) 
behavior: In the first one [||,||J|] it is assumed that dis- 
order introduces a distribution of (one-channel) Kondo 
temperatures Tr in the system; this distribution can be 
directly related to an anomalous low temperature behav- 
ior. The second model proposes a single-impurity origin 
of the NFL, e.g. associated with the quadrupolar (two- 
channel) Kondo effect 
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FIG. 1. Phase diagram of CeCue-aiAu^ as a function of 
doping and temperature. The points are Neel temperatures 
determined by different methods |7j] (open and closed sym- 
bols for single and polycrystals respectively), the solid line 
denotes the phase transition, the dashed lines are theoretical 
crossover lines. The regions are described in the text. 

In CeCu6- K Au K , there is clear experimental evidence 
[^-0 for a third mechanism based on the proximity to a 



quantum phase transition (QPT) [p~0| jl3|] near x c = 0.1. 
For x > x c the system orders magnetically with a Neel 
temperature TV oc (x — x c ) M with jj, = 1 as shown in 
Fig. |l|. The QPT can be interpreted as the result of 
the competition between the Kondo effect, which tends 
to screen the magnetic moments, and the RKKY inter- 
action, which favors a magnetically ordered state. The 
Kondo effect is weakened by increasing doping which, as 
suggested experimentally, only leads to a volume increase 
with no apparent change in the number of carriers. 
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FIG. 2. The specific heat c/T of CeCue-x Au^ versus log T 
(from @). 

CeCug-zAuz remains metallic for all x with typical FL 
properties at low temperatures both at x = and for x ^> 
x c . However, at x — x c , Tn vanishes and NFL behavior 
is observed in all accessible quantities down to the lowest 
temperatures, T <C Tk ~ 6if . As an example we show in 
Fig. the specific heat with c/T oc ln(To/T) over nearly 
two decades. At the same doping, the static susceptibility 
X oc 1 — aVT and the resistivity p ~ pa + A'T show a 
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remarkable NFL behavior over a substantial T range || . 

Up to now [|U^,|lJ] it was assumed that the critical 
fluctuations of the QPT are dominated by incommensu- 
rate three-dimensional correlations. This would strongly 
suggest a description by a quantum-critical theory with 

LI 14 1 . However, 



d = 3, z = 2 as investigated in ||K 
this well-established theory is in contradiction to the 
thermodynamic properties which are observed in the ex- 
periment, as it would suggest c/T oc 1 - ByT and 
p p + A"T 3/2 at low temperatures One could, 

however, argue that in the experiments only a crossover 
region is accessible JL4|. In our opinion this is not fully 
convincing and it appears not to be possible to fit both 
resistivity and specific heat over the observed T range by 
such a theory 0,^5) . 
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FIG. 3. Elastic neutron scattering of CeCu5.8Auo.2 along 
the a axis at 70mK. The inset shows that all features are 
sharp in the b direction. 

Below we will focus on a novel feature of the magnetic 
fluctuations which we uncovered by carrying out elastic 
neutron scattering experiments near the quantum criti- 
cal point (QCP). The data taken in the ordered phase 
(x = 0.2) is displayed in Fig. [| A scan along (ft, 0, 0) 
reveals besides the nuclear reflections two magnetic fea- 
tures which vanish for higher temperatures. The satel- 
lite peaks at (±0.79, 0, 0) describe a three-dimensional 
incommensurate magnetic order. The magnitude of the 
ordered magnetic moment is extremely small and is es- 
timated to be of the order of 0.02/is. The correlation 
length determined from the peak width is approximately 
20 lattice constants a (orthorhombic notation), with 4 Ce 
atoms per unit cell. 

The second feature, seen as a broad sine-modulated 
background with maxima at (1, 0, 0) and (3, 0, 0) in Fig. 
||, has a factor of 3 — 4 higher integrated weight. It is 
important to stress that both features are sharp in the 
b direction with a width of order 20 unit cells as can be 
seen from scans along (ho,k, 0) for different ha (Fig. ||). 



In this letter we propose that the magnetic fluctuations 
associated with the second broad structure in the neutron 
scattering dominate the critical fluctuations at the QPT 
in the observed temperature range. Starting from this 
assumption we will show that the logarithmic increase of 
c/T for decreasing temperature, the linear resistivity and 
the phase diagram can easily be derived. We interpret 
this structure as arising from ferromagnetically ordered 
planes perpendicular to the a direction. The sine mod- 
ulation would suggest that effectively two of these ferro- 
magnetic planes with distance a/2 couple antiferromag- 
netically in the a direction, but different pairs of planes 
are incoherent. While it is possible to identify slightly 
corrugated planes of Ce atoms in the crystal structure, 
we do not see an obvious reason for a strongly asymmetric 
coupling, which could directly explain the observed two- 
dimensional structure. For the following, the details of 
the magnetic order are irrelevant and we will only assume 
that some critical two-dimensional fluctuations exist. 

To derive an effective action near the phase transition 
in a Ginzburg-Landau- Wilson approach we have to 
account for the damping of the critical fluctuations. To 
describe the critical modes, we introduce a scalar field 
$ q . • the experiments suggest a preferred direction of the 
magnetic moments along the c axis [|7|]l6|]. q|| is a two- 
dimensional vector in the be plane. Ordering will occur 
at q|| = 0. The primary damping mechanism is the cou- 
pling to particle-hole pairs. The dynamics of the quasi- 
particles is three-dimensional, this can be inferred from 
the transport properties which vary at most by a factor of 
two in different directions jl5| . We assume a coupling of 
the critical fluctuations to the heavy quasi-particles (with 
creation operators c^) by the following Hamiltonian: 
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(1) 



q = q|| + q± is the transferred momentum split up 
in two components parallel and perpendicular to the 
planes, a 1 are the Pauli matrices and g is the cou- 
pling constant. h(q±) is some smooth function describ- 
ing the magnetic structure perpendicular to the planes, 
e.g. h(q±) = isinqj_a/4 for two planes with distance a/2 
which are coupled antiferromagnetically. The details of 
H c are, however, not important for our discussion. 

Integrating out the fermions (cf. [lC|]) induces damping 
of and interaction between the critical modes ^q^,^- To 
quadratic order, we obtain the following contribution to 
the effective action in imaginary time (/3 = \/ {ksT)): 



Wn,q|| 



Q| • 



uj n = 2nn/P are bosonic Matsubara frequencies. The 
damping is described by the imaginary part of the 
particle-hole bubble x°> I™ J dq±\h(q±)\ 2 \° {q, w + «0) ~ 
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ju>. Near the QCP, after a proper rescaling, the effective 
action takes the following form: 



S — S*2 + Sint 
1 x - 



(2) 



S 2 = 
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iS'int describes the interaction between the critical modes, 
the leading term is given by S nlt « U J^dr J G? 2 r|$(r, r)| 4 
p0| . The distance from the QCP is measured by S. At a 
critical value S Cl determined by the strength of interac- 
tion, the system is at the QCP, i.e. S — 8 C oc x c — x. The 
effective action corresponds to a quantum-critical theory 
in d = 2 with a dynamic exponent z = 2. This dynamic 
exponent describes the fact that if one scales momenta 
by k — * Ak one has to scale frequencies by u — ► X z ui. 

This theory has been widely investigated by many au- 
thors (e.g. []l0|-|l2"|| ) , mainly in the context of the anti- 
ferromagnetic spin -fluctuation picture of high-T c com- 
pounds. We will primarily employ the results of Millis 
Jll| who did a careful renormalization-group study of Eq. 
(0). The effective dimension of this theory at T = is 
d + z = 4; the scaling dimension of the interaction Si n t 
is 4 — (d + z) = 0, therefore it is marginal. Nevertheless, 
the leading behavior of the specific heat in the disordered 
phase can be directly calculated from the Gaussian part, 
52, of the action of Eq. (0). The free energy per volume 
corresponding to S2 is 
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Note that <J C = in this approximation; oj c is a typical 
cutoff energy of the order of the Kondo temperature in 
the system. Consequently, the coefficient of the specific 
heat 7 = c/T = —d 2 F/dT 2 diverges logarithmically with 
decreasing temperature at the QCP, as observed in the 
experiment. For 5 > 5 C , i.e. x < x c , 7 stays finite. 
Note that the logarithmic terms do not arise from some 
marginal operators, they are due to a pure phase space 
effect, this typically happens at d + z — 2z EJ. 

By a solution (Tl]] of the scaling equations for the model 
of Eq. (0) a phase diagram emerges as displayed in Fig. 
[l]. Region I, the low temperature phase for x > x c , is 
the ordered phase. The behavior near Tn will depend on 
the structure of the order parameter, e.g., whether the 
ordering is of Ising type or how it is stabilized by three- 
dimensional coherence. The Ginzburg criterion for our 
model in Eq. (0) predicts |ll[] that Tn is proportional 
to the distance from the critical point, Tn oc x — x c . 
This relation is fulfilled by the experiment over an aston- 
ishingly large range. The other crossover lines are hard 
to analyze quantitatively with the existing experimen- 
tal data, therefore we do not attempt to fit them. The 



qualitative trend is, however, consistent with the phase 
diagram. In Fig. |f| we include the theoretical curves 
of Q. Re gion 11 is dominated by the fluctuations of 
the finite-temperature phase transition; in region III true 
quantum-critical behavior with 7 oc ln(To/T) and, as we 
shall see, p ps po+A'T can be observed. In this quantum- 
critical regime - in [^f) it is called "classical Gaussian 
regime" - the temperature is the most important energy 
scale. The correlation length £ depends on temperature 
as oc T with logarithmic corrections. In region IV, a 
pure crossover regime which is probably hard to observe, 
£ is determined by the distance from the critical point, 
oc x c — x, while energy fluctuations and therefore 
the specific heat are governed by the temperature with 
7 oc ln(To/T). For region V, FL behavior is expected 
with a linear specific heat c/T oc — ln(x c — x) and a finite 
correlation length £~ 2 oc x c — x as in region IV. 

The whole scenario is in qualitative agreement with the 
experiments In particular, the logarithmic increase 

of the specific-heat coefficient in region III, the T-linear 
resistivity which we will calculate in the following, and 
the linear increase of Tn with x are confirmed by the 
experiments. 

For the calculation of the resistivity we closely follow 
Hlubina and Rice jl7| who have calculated the resistiv- 
ity of electrons in two dimensions coupled to quantum- 
critical antiferromagnetic spin fluctuations. In our case 
the dynamics of the quasi-particles is three-dimensional 
as mentioned above. Our starting point is Eq. (|f|) and 
the corresponding collision term in a Boltzmann equa- 
tion: 
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Here we have already linearized the collision term; the 
occupation of a state with momentum k is given by 
/k = /£ + ^(df^/de) where /° = /°(e k ) is the usual 
Fermi function; n(e) is the Bose function. We have 
omitted the factor \h(q±)\ 2 , the qualitative behavior of 
the resistivity is not influenced by this smooth function. 
Xq( w ) = u$i,ui) is the order-parameter susceptibil- 
ity. In the following we use for the disordered phase 
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T* + cT + q 



(5) 
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with temperature-independent constants A,c and T*. 
This phenomenological form corresponds to the behav- 
ior of the correlation length described above, £~ 2 oc T 
in the quantum-critical region up to logarithmic correc- 
tions. T* oc S — 5 C measures the distance from the QCP. 
Eq. (||) is valid in regions III-V, we are primarily inter- 
ested in the quantum-critical regime III with T* <C T. 
Following fl7||, the resistivity can be determined within 
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Boltzmann theory from the minimization of a functional 
of v? k 
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where h is the direction of the applied electric field, v k 
the velocity of the electrons, po — h/e 2 and 



/k'M e k' - £k)ImXk' -kn («*' ~ £k). 



Note that a second contribution due to impurity scatter- 
ing has to be added, which is not given here. For sim- 
plicity we assume a spherical Fermi surface. The radial 
part of the momentum integration can be rewritten as 
an energy integration. Using J /°(e)(l — /°(e + ui))de = 
lo(1 + n(u>)) and f Q um(u))[l + n(w)]Imxk> -k N (u)dw = 
I[f3{T* + cT + (k^ - k||) 2 )]TA with I[x] « tt 2 /(3x(x + 
2n/3)) [0 we can perform the energy integration. As 
long as the resistivity is dominated by impurity scatter- 
ing - this is true in the whole range where a linear tem- 
perature dependence has been observed - we can assume 
p7[ <y9 k ~ v k • n and arrive at 



ApcxT 2 
7(«) = 
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(T* + cT + k 2 F a)(T* + c'T + k 2 F a 
dfi k dQ k >(n ■ k - n- k') 2 S(a - (kf, - k,,) 2 ) 



(||) gives l/r k cx e k In 1 /e k for directions parallel to the 
planes and l/r k cx e k in the o direction. Averaging over 
the Fermi surface results in (l/rk) oc e k lnl/e k , in sharp 
contrast to the usual l/r k cx e 2 . in a Fermi liquid. 

In |&R the static susceptibility for x = x c was fitted 
by % cx 1 — ol\[T from the lowest measuring tempera- 
ture of 80 mK up to 3 K. However, the data can also 
be well described by x ~ °o + l/(ai + CL2T) for temper- 
atures up to 1.4 K. One would expect a susceptibility 
of the last form for two antiferromagnetically coupled 
planes with an order-parameter susceptibility as given 
by Eq. (@) . Further theoretical and experimental studies 
are needed to clarify this point. It will also be important 
to investigate further the interplay of the two- and three- 
dimensional order which could finally lead to a change 
in the observed properties at some lower temperature. 
We think that the explanation of the phase diagram, of 
the linear increase of the resistivity with temperature 
and of the anomalous specific heat in CeCu6-2;Au x by 
a quantum-phase transition with d = 2, z = 2 is already 
promising. To our knowledge, this would establish the 
first clear experimental realization of such a theoretical 
scenario. 
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ported by the DFG and LIP. 



with c' — c + 2tt/3. k denotes a unit vector in the direc- 
tion of k and k|| its projection to the be plane. For n not 
parallel to the planes, we find 17(a) ~ const for small a 
and therefore we get approximately 



T , T* + cT 
Ap cx In 



T* + c'T 



T 2 ,T<T*/c' 

In c/c' , . (8) 

T '— ,T> T*/c y ' 



For a finite T* and small temperatures we recover the 
usual Ap cx T 2 of a FL. However, in the quantum crit- 
ical regime, i.e. for T* <C cT, the resistivity is linear in 
temperature as observed in the experiment (^J^,^| . 

For an electric field parallel to the planes, we obtain 
g(a) cx —a In a and accordingly Ap cx T 2 lnT for T* = 0. 
This is however an artifact of our approximation of using 
a spherical Fermi surface and an isotropic s-wave scat- 
tering amplitude. In a more realistic approach, different 
directions would mix and give a linear increase of re- 
sistivity with temperature in all directions, as observed 
M|. Nevertheless, one would still expect that the lin- 
ear increase of the resistivity at the QCP is largest in 
the direction perpendicular to the plane which is indeed 
observed p"5[ . 

The NFL character is also manifest in an anomalous 
self-energy of the electrons. Calculation of the lifetime 
in Born approximation at T — using Eqs. (0) and 
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